Functional superconductor interfaces from broken time-reversal symmetry 
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The breaking of time-reversal symmetry in a triplet superconductor Josephson junction is shown 
to cause a magnetic instability of the tunneling barrier. Using a Ginzburg-Landau analysis of the free 
energy, we predict that this novel functional behaviour reflects the formation of an exotic Josephson 
state, distinguished by the existence of fractional flux quanta at the barrier. The crucial role of 
the orbital pairing state is demonstrated by studying complementary microscopic models of the 
junction. Signatures of the magnetic instability are found in the critical current of the junction. 
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The Josephson effect between superconductors sepa- 
rated by a tunneling barrier continues to be of fundamen- 
tal interest, in particular as a phase-sensitive test of the 
pairing symmetry of unconventional superconductors)^ 
and in the study of the interplay of superconductivity 
with magnetism.™!! Although the qualitative features of 
the Josephson effect are determined by the quantum na- 
ture of the superconductors, the modification of the su- 
perconducting state at the junction interface must often 
be included in the proper description of the supercurrent 
transmission, the so-called proximity effect lii^i! In con- 
trast, the properties of the barrier are usually regarded 
as fixed.— Recently, however, it has been shown that a 
thin ferromagnetic layer on a singlet superconductor can 
display novel behaviour."* In particular, the presence of 
the superconductor can suppress the uniform magnetiza- 
tion^ or stabilize a domain structure in the ferromagnet?^ 
These effects result from the competition between singlet 
superconductivity and ferromagnetism, and indicate that 
the tunneling barrier in a Josephson junction is not nec- 
essarily independent of the superconductors. 

In this letter, we consider the possibility of using the 
presence of the two superconductors to induce a mag- 
netic instability of a non-magnetic tunneling barrier in 
a Josephson junction. In particular, by both general 
phenomenological arguments and solution of specific mi- 
croscopic models, we show that such a novel functional- 
ity of the barrier can develop for time-reversal symmetry 
(TRS) breaking configurations of two spin-triplet super- 
conductors on either side. The Josephson coupling across 
the tunneling barrier is essential to this effect, which 
manifests itself as an exotic state distinguished, for ex- 
ample, by the existence of fractional flux quanta at the 
barrier.^ Moreover, such a junction displays an anoma- 
lous temperature dependence of the critical current. 

The intrinsic spin structure of the Cooper pairs in a 
triplet superconductor (TSC) requires a quasiparticle gap 
with three components Ag^ for each z-component of spin 
Sz ~ —1,0, +1. Each of these gaps has odd orbital par- 
ity, i.e. A5^(— k) = — A5^(k). The pairing state is con- 
veniently described by the so-called d-vector, defined in 



spin-space d = ^(A^i - Ai)x - |(A_i + Ai)y + Aqz, 
which also serves as the order parameter for the TSC. 
Although a multitude of different triplet pairing states 
are allowed by symmetry, only a few examples of TSCs 
have been discovered so far, e.g. Sr2Ru04j^ UGe2i^ In 
Sr2Ru04 the spin pairing state has been identified as 
unitary and equal-spin-pairing,*- i.e. the spins of the 
triplet Cooper pairs lie in the same plane in spin space 
perpendicular to d, but the condensate does not have a 
net spin. Restricting ourselves to such pairing states, we 
write d — de**^ where d is a real vector. The vector char- 
acter of the TSC order parameter provides a novel degree 
of freedom in Josephson junction physics: in addition to 
the phase difference between the condensates to the left 
and right of the barrier, which controls the Josephson 
supercurrent as in the familiar spin-singlet case, there is 
also the mutual alignment of the left (d^) and right (d/?) 
vectors, which controls the magnetic aspects of the trans- 
port.** In particular, when d^, x d/j 7^ 0, a Cooper pair 
tunneling across the barrier undergoes a reconstruction 
of its spin state, producing an effective spin through the 
TRS breaking combination (S) = id^ x d^ + H.c. at the 
junction interface. 

The violation of TRS at the tunneling barrier by the 
misaligned d-vectors allows the TSCs to directly couple 
to a magnetization M of the interface. As we will see, 
the TSCs may in fact change the electronic properties of 
the interface so as to stabilize a spontaneous ferromag- 
netic order. This novel behavior can be understood on 
a phenomenological level by a Ginzburg-Landau analysis 
of the free energy. Introducing the TSC order parame- 
ters for each side of the interface as d^ = d^e^'^^ and 
dji = dflc*"^", we write the free energy of the junction F 
to lowest order in M and d^ jj as 



F 



IMP ~ ~ ~ ~ 

^ - tdi • dR cos{(b) + 27M • (di X dfl) sin(,^) (1) 



Here (p — <j)R — (j)L is the phase difference, x denotes the 
intrinsic uniform spin susceptibility of the barrier, and 
t and 7 are phenomenological parameters. The ground 
state of the junction is obtained by minimizing F with 



triplet ferromagnetic triplet 

superconductor barrier superconductor 




z<0 



!>0 



FIG. 1: (color online) (a) The Pz-Pz junction, with pz-wave 
orbital pairing in the two superconductors, d-vectors in the 
x-y plane, and an induced magnetic moment along the z- 
axis. (b) The Py-py junction, with py-wave orbital pairing in 
the two superconductors, d-vectors in the x-z plane, and an 
induced magnetic moment along the y-axis. 



respect to M and (f>. We find that the non-magnetic state 
of the barrier is unstable for xi'^l^L x dj^p > Ml ■ dj^: 
the coupling to the intrinsic spin (S) = 2dj;, xd/j sin(<^) of 
the junction instead realizes a TRS breaking state, char- 
acterized by a non- vanishing magnetization M || d^ x d^^, 
and a phase 4> = 4>min 7^ 0, tt. The junction is then in an 
exotic fractional stated where the magnetic barrier is ca- 
pable of carrying flux lines with non-integer multiples of 
the flux quantum $0 = hc/2e. The observation of this 
characteristic feature is discussed below. 

The magnetic instability depends essentially upon, and 
can be tuned by, the mutual misalignment of the d- 
vectors of the two TSCs. It is also controlled by the 
specific details of the junction, through the susceptibility 
X, and the parameters t and 7. The latter are fixed by 
the orbital part of the bulk pairing state. To elucidate 
the crucial role this plays in the magnetic instability, we 
examine two complementary microscopic models of the 
junction. The first, shown in Fig. ©(a), has the TSCs in 
a p^-wave orbital state (the Pz-Pz junction), and the left 
and right d-vectors are aligned parallel to the barrier but 
at an angle 2ri with respect to each other. In the second 
model the orbital state is py-wave (the Py-Py junction), 

and the d-vectors lie in the x-z plane but are again mis- 
oriented by the angle 2ri [Fig. (II|)(b)]. We assume trans- 
lational invariance in the x-y plane and that the bulk 
TSCs extend indefinitely along the z-axis. Furthermore, 
we take the TSCs to have spatially-constant gaps, and 
the maximum gap magnitude A displays weak-coupling 
temperature-dependence, with T = value Aq. The tun- 
neling barrier is modeled to be of (5-function width, with 
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FIG. 2: (color online) The Andreev bound state spectrum in 



the (a) Pz-pz and (b) the Py-Py junction for kz 



ky 



kpl^Pl. 



normal-state height Up = Zhvp where Z is a dimen- 
sionless quantity and vf = hkF/rn is the Fermi velocity 
in the bulk TSCs, which are assumed to have spheri- 
cal Fermi sufaces of radius kp- If the barrier supports 
a magnetic moment M || d^ x dji, the effective barrier 
height for spin-cr quasiparticles with spin parallel to M in 
dimensionless units is Z — aM where M = gfisl^l/fivp. 

The electronic properties of the two junctions can be 
expressed entirely in terms of the Andreev bound state 
(ABS) spectrum.^ These subgap states are localized at 
the barrier interface and are formed by multiple Andreev 
refiection of tunneling quasiparticles within the barrier. 
We solve the Bogoliubov-de Gennes equations^- to obtain 
explicit expressions for the ABS energies Ei^^^ in the Pz- 
Pz junction 



Sk,. = \A{T)kz\V%{k) cos(</)/2 - ai]) 



and the Py-Py junction 



(2) 



Eu^, = \A{T)ky\^l-%{k)sin'i(f>/2 + aT]) (3) 

The bound states are indexed by the component of spin 
(T = ±1 parallel to d^ x d^, and each state has two 
branches at ±E^. The transparency of the barrier to 
spin-(T quasiparticles is given by 7^(1^) = fc^/[fc^ -\- (Z — 
aM)'^kp]. We plot the ABS spectrum as a function of (j) 
in Fig. (21). Note that the bound states always intersect 
the line E — in the Pz-Pz junction. These so-called zero 
energy state are guaranteed by the arrangement of the p- 
wave orbitals, such that all specularly-refiected quasipar- 
ticles experience a sign-change of the gap. Since the gap 
does not change sign for refiected quasiparticles in the 
Py-Py junction, in contrast, zero energy states are only 
found here for a perfectly transparent tunneling barrier, 
as is also the case for s-wave superconductor junctions^ 
The electronic contribution to the free energy of the 
junctions can be written in terms of the ABS energies 



2 cosh 



2kBTj 



(4) 



As in Eq. ([T]), the magnetic free energy of the barrier 

M^/2x, where x is 



is included to lowest-order F, 



mag 



given in units of {g^s/f^vp) / Aq. Numerically minimiz- 
ing the total free energy F^i -\- Fmag with respect to both 
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FIG. 3: (color online) Induced magnetic moment M and sta- 
ble phase difference c/jmin as a function of reduced temperature 
for (a) the Pz-Pz and (b) the Py-Py junction. In both panels 
we set ?7 = 0.2n, Z = 0.7 and x ~ 20. Magnetic phase dia- 
gram for the (c) Pz-Pz and (d) Py-Py junctions as a function of 
Z and T. A magnetic moment is stable in the region labeled 
FM, while the region of non-magnetic behaviour is denoted as 
PM; metastable states are shown in brackets. Second-order 
transitions are indicated by a solid line, first-order transitions 
by a dashed line, and the limits of the metastable states by a 
dotted line, rj and x are as in (a) and (b). 



M and 0, we find the global free energy minimum. Typ- 
ical minimizing values for the Pz-Pz and Py-Py junctions 
are plotted as a function of temperature in Fig. ©(a) 
and Fig. ©(b) respectively. We find that the barrier 
undergoes a magnetic instability at sufficiently large Xj 
and below a critical temperature Tm < Tc such that the 
magnetic state appears only in the presence of super- 
conductivity. For Z ^ the junction is in a fractional 
state below Tm with two degenerate free energy minima 
(M,0„iin) and (-A/,-0mi„) (broken TRS). 

In Fig. ©(c) and Fig. (I3])(d) we show the phase dia- 
gram as a function of Z and T at fixed x for the Pz-Pz and 
Py-Py junctions respectively. The qualitatively different 
form of these phase diagrams follows from the response 
of the ABS spectrum to the appearance of the magne- 
tization and the resulting change in the 7^(k). As can 
be seen from Eq. ^ and Eq. I©, the ABS spectrum 
in the two junctions has very different dependence upon 
7^(k): due to the zero-energy states in the Pz-Pz junc- 
tion, each bound state E^u monotonically shifts towards 
the middle of the gap with decreasing 7^(k); for the Py- 
Py junction, in contrast, the states move towards the gap 
edges. From Eq. (|4]), the free energy contributed by £^k,(T 
in a Pz-Pz junction will therefore increase (decrease) as 
the transparency 7^(k) decreases (increases), while the 
opposite is true for the Py-Py junction. 

It immediately follows that the Pz-Pz junction is non- 
magnetic at Z = 0, as M 7^ would reduce the trans- 
parency for both spin orientations and hence raise the 




FIG. 4: (color online) Proposed experiment for the observa- 
tion of fractional fiux quanta at the interface between mag- 
netic and non-magnetic regions of a tunneling barrier. The 
contour C is used in evaluating Eq. ([Sjl. 



total free energy. The M = state remains stable be- 
low some critical value of Z; increasing Z beyond this, a 
magnetic moment appears as the free energy gain from 
decreasing the transparency in the a = —1 sector out- 
weighs the increase in the <t = +1 sector. Although 
the decrease in electronic free energy favors the indefinite 
growth of M with increasing Z, the maximum magnitude 
of M is limited by the cost in magnetic free energy. For 
sufficiently large x and low temperatures, the transition 
back into the non-magnetic state is first order, with re- 
gions in the phase diagram where the non-magnetic and 
magnetic states are metastable, as shown in Fig. ([31) (c). 
At higher temperatures or smaller Xi \M\ continuously 
vanishes after going through a maximum. 

In contrast, the Py-Py junction displays a spontaneous 
magnetization at Z — for all T < T^, stabilized due to 
the reduction in Fei from the decreased transparency in 
each spin sector. This effect is absent from the Ginzburg- 
Landau expansion of F^i in Eq. ([T]), as we have only kept 
terms to first order in M; from Eq. ([3]), however, we 
see that the magnetization only enters Fei as |Mp when 
Z = 0. Despite the spontaneous magnetization, the junc- 
tion is not in a fractional state and </'min = 0. Turning 
on a finite tunneling barrier strength (Z > 0) at fixed 
T, the magnetic moment of the barrier decreases to com- 
pensate for the free-energy increase from the enhanced 
transparency in the a = +1 sector; 4'rain simultaneously 
takes on a fractional value. As Z is further increased, 
the barrier moment is monotonically suppressed, while 
the stable phase difference passes through a stationary 
point before returning to its Z = value. 

The characteristic signature of the magnetic instability 
is the appearance of fractional flux quanta at the junc- 
tion interface. In Fig. ^ we show a proposal for their 
observation in a Josephson junction with a tunneling bar- 
rier consisting of two materials, one of which undergoes 
the magnetic instability proposed here, while the other 
remains non-magnetic at all temperatures. The stable 
phase difference across the magnetic region is 0, while 
it is across the non-magnetic region. If a magnetic 
flux line is trapped at the boundary between the barrier 
segments, a line integral along the contour C in Fig. Q 
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FIG. 5: (color online) (a) Current vs phase relationships in the 
Py-Py junction at T = 0.2Tc both with and without (Af — 0) 
the magnetic instability, (b) Critical Josephson current as 
a function of reduced temperature in the Py-Py junction. In 
both panels we take rj = 0.2tt and Z = 0.7. 



shows that the enclosed flux $ is 



— = n+ (i)ds-V0 = n+ — , n eZ (5) 
$0 Jc 27r 

As the temperature is lowered below Tm, 4> takes a frac- 
tional value, and a flux line appears with a continuously 
increasing flux <I> 7^ $o- Experimentally, such a flux line 
could be directly observed by local magnetic probes like 
scanning SQUID microscopy, or inferred from the asym- 
metric Fraunhofer pattern of critical current vs applied 
field. This proposal resembles the devices used to ob- 
serve half- integer flux quanta by Weides et al.^ where 
spin-singlet superconductors were used instead of spin- 
triplet, and differing widths of a permanent magnetic 
barrier guaranteed (p — n always. Although other pro- 
posals exist for the creation of fractional flux quantaj^ii^ 
their detection in our proposed junction would be unam- 
biguous confirmation of the magnetic instability. 

The magnetic instability of the tunneling barrier rad- 
ically alters the supercurrent transmission through the 
junction. From Eq. ([l]), we find that the Joseph- 
son current vs phase relationship // = {e/fi)dF/d<j) oc 
sin((/) — (/)niin) is shifted from its usual form for a non- 



* Electronic address: |brydon@theory.phy.tu-dresden.de| 

^ S. Kashiwaya and Y. Tanaka, Rep. Prog. Phys. 63 1641 
(2000). 

^ C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

^ A. A. Golubov, M. Y. Kupriyanov, and E. Il'ichev, Rev. 
Mod. Phys. 76, 411 (2004). 

* A. I. Buzdin, Rev. Mod. Phys. 77, 935, (2005); F. S. 
Begeret, A. F. Volkov, and K. B. Efetov, ibid 77, 1321 
(2005). 

^ F. S. Begeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 

B 69 174504 (2004). 
^ A. I. Buzdin and L. N. Bulaevskii, Sov. Phys. JETP 67, 

576 (1988); F. S. Bergeret, K. B. Efetov, and A. I. Larkin, 

Phys. Rev. B 62, 11872 (2000). 



magnetic barrier. This is clearly seen in Fig. ([5]) (a) for 
the Py-Py junction (the Pz-Pz junction results are qual- 
itatively identical). Note that the different magnitudes 
of max{/j} and min{/j} are due to higher-order har- 
monics in (/) which are not included in the free energy 
expansion Eq. ([T]). We also find a strong enhancement 
of the critical current if — max{|/j|} below the mag- 
netic instability [Fig. (in])(b)]. This occurs as the in- 
creased current through the spin sector with the en- 
hanced transparency over-compensates for the decreased 
current through the spin sector with the lowered trans- 
parency. The increase of below Tm is reminiscent 
of the "low-temperature-anomaly" of d-wave Josephson 
j unctions fi 

For simplicity, we have neglected the suppression of the 
TSC state near the interface due to the proximity effect. 
Including this would increase F^i and hence shrink the 
parameter space where the moment is stable. Although 
quantitative changes in the phase diagrams Fig. (l3|)(c) 
and (d) are expected, the free energy expansion Eq. ([T]) 
and the basic experimental signatures of the fractional 
state remain valid. Furthermore, we do not expect a 
qualitatively different ABS spectrum, and so the essential 
role of the orbital pairing state of the two TSCs should 
remain in a fully self-consistent analysis. 

In conclusion, we have shown that a TRS-breaking con- 
figuration of two TSCs in a Josephson junction can cause 
the tunneling barrier to develop a spontaneous magne- 
tization. This realizes an exotic Josephson state with 
stable phase difference < 0niin < t^- The orbital part 
of the TSC pairing state was demonstrated to control 
the magnetic instability. The existence of fractional flux 
quanta at the barrier, and a large increase in the criti- 
cal current beneath the magnetic transition temperature, 
are the experimental signatures of this state. 
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